System and method for medical treatment hypothesis testing

ABSTRACT

A computer-implemented system and method of evaluating the effects of medical treatments, the method including receiving patient record data; identifying relevant characteristics for evaluation; identifying a first treatment; identifying a second treatment; assigning a weight to each patient case; determining the relative likelihood, using the assigned weights, that an identified treatment will result in an identified effect when contrasted with a second identified treatment; and, outputting this estimated relative likelihood.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application is a continuation application of U.S. application Ser. No. 12/815,356, filed Jun. 14, 2010, which claims the benefit of U.S. Provisional Application No. 61/186,759, filed Jun. 12, 2009, the contents of which are expressly incorporated herein by reference in its entirety.

FIELD OF THE INVENTION

The present invention relates generally to methods and systems for the statistical analysis of retrospective medical data and more specifically to methods and systems for the evaluation of the effects of medical treatments.

BACKGROUND OF THE INVENTION

There exist systems and methods for analyzing the effects of medical treatments using cohort studies. Most of these systems and methods are based on randomized controlled trials of particular treatments. Given a sufficient number of subjects over a sufficient period of time, randomized controlled trials have the advantage of simplifying this analysis by evenly distributing confounding factors (or differences between the groups) across the group receiving the treatment and the control group. These systems and methods are limited because randomized controlled trials are expensive to conduct, they rely for their accuracy in part on a high number of test subjects, and they are not very effective at identifying secondary interactions resulting from conditions that are a-typical among treatment candidates. Additionally, randomized control trials, which are prospective cohort studies, cannot capitalize on the wealth of information available in the broader set of existing patient record data in the same way that retrospective cohort studies can.

Other existing systems and methods analyze historical patient record data to evaluate the effects of particular treatments. These systems are also limited because they use a largely manual modeling process to control for confounders and patient characteristics that are not evenly distributed in the historical data across the treated group and non-treated groups. Because of the high risk of bias in these manual confounder control processes, the quality of the results produced by the existing systems are highly dependent on the operator's expertise level. Additionally, main and interaction effects associated with a high number of potential confounders in large existing data sets makes manual modeling very time intensive and subject to human error.

SUMMARY OF THE INVENTION

The invention relates to a software component system and method for the statistical analysis of patient medical records for the purpose of evaluating the relative likelihood that patient experience identified effects when treated using identified treatments. This method of analysis produces results that decision makers in the medical field can rely on to evaluate the efficacy and risks of identified treatments in determining which treatments to approve, use, and fund.

In some embodiments the system evaluates two treatments to determine the relative likelihood that patients receiving the first treatment will experience the identified effect by comparison to patients receiving the second treatment. In other embodiments, the treatments being compared are actually dosage amounts of the same treatment. In yet other embodiments, the system compares a first treatment to several different second treatments by performing several individual comparisons to produce several different ratio of rate ratios, each representing the relative likelihood that the patient will experience the identified effect as between the first treatment and another treatment.

One embodiment of the present invention provides a computer-implemented system and method of evaluating the effects of medical treatments, the method including receiving patient record data, identifying relevant characteristics for evaluation, identifying a first treatment, identifying a second treatment; assigning a weight to each patient case, determining the relative likelihood, using the assigned weights, that an identified treatment will result in an identified effect when contrasted with a second identified treatment, and outputting this estimated relative likelihood.

In another embodiment the invention relates to a system for evaluating the effects of medical treatments including: a network interface, a patient record database residing on a server accessed through a network, and a server for analysis, wherein the analysis server is configured to receive patient record data, identify relevant characteristics for evaluation, identify a first treatment, identify a second treatment, assign a weight to each patient case based on the likelihood that the patent would be a member of the exposed group, determine the relative likelihood, using the assigned weights, that an identified treatment will result in an identified effect when contrasted with a second identified treatment, and, output this estimated relative likelihood.

In another embodiment the invention relates to a server for evaluating the effects of medical treatments, the server including: a processor; and memory operably coupled to the processor storing programming instructions therein, the processor being operable to execute program instructions, the program instructions including: assigning a weight to each patient case based on the likelihood that the patent would be a member of the exposed group; using the assigned weights in determining the relative likelihood that an identified treatment will result in an identified effect when contrasted with a second identified treatment.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a schematic block diagram of a medical treatment hypothesis testing system constructed according to an embodiment of the present invention for comparing the effects of medical treatments based on retrospective observational data;

FIG. 2 is a high-level flow diagram of the steps exercised by the data selection component of the medical treatment hypothesis testing system of FIG. 1 according to an embodiment of the present invention;

FIG. 3 is a high-level flow diagram of the steps exercised by the propensity scoring component of the medical treatment hypothesis testing system of FIG. 1 according to an embodiment of the present invention;

FIG. 4 is a flow diagram of the steps exercised by the Prepare Data in Table component of the propensity scoring component of FIG. 3 according to an embodiment of the present invention;

FIG. 5A is an illustration of an exemplary patient case data table at the beginning of the Prepare Data in Table process of FIG. 4 according to an embodiment of the present invention;

FIG. 5B is an illustration of an exemplary patient case data table at one point in the Prepare Data in Table process of FIG. 4 according to an embodiment of the present invention;

FIG. 5C is an illustration of an exemplary patient case data table at one point in the Prepare Data in Table process of FIG. 4 according to an embodiment of the present invention;

FIG. 5D is an illustration of an exemplary patient case data table at one point in the Prepare Data in Table process of FIG. 4 according to an embodiment of the present invention;

FIG. 5E is an illustration of an exemplary patient case data table at one point in the Prepare Data in Table process of FIG. 4 according to an embodiment of the present invention;

FIG. 6 is an illustration of an exemplary table for the comparison of non-weighted patient data for patients in the treatment group and weighted patient data for patients not in the treatment group according to an embodiment of the present invention;

FIG. 7 is an illustration of an exemplary patient case data table for use in the regression modeling 118 of FIG. 1. This illustration shows the data as it is formatted during the format data for regression step 312 of FIG. 3.

The figures depict embodiments of the present invention for purposes of illustration only. One skilled in the art will readily recognize from the following discussion that alternative embodiments of the structures and methods illustrated herein may be employed without departing from the principles of the invention described herein.

DETAILED DESCRIPTION OF THE INVENTION

Turning now to the drawing, embodiments of the present invention that include patient data selection, propensity scoring, inverse-probability weighting, and doubly robust estimation of the ratio of rate ratios are shown. Retrospective observational medical record data are selected by an operator from a patient record database based on two identified treatments and identified effects to be studied. Each patient who did not receive the first identified treatment is assigned a propensity score that represents the likelihood, based on the patient's individual characteristics, that the patient would have been in the group that received the first identified treatment. The assigned propensity scores are used to weight each patient's data so that the weighted data of the group that did not receive the first identified treatment (i.e. the group that received the second identified treatment) closely resembles the un-weighted data of the group that received the first treatment. The weighted data is used to perform a regression to estimate the relative likelihood that a patient receiving the first identified treatment would experience an identified effect as compared to a patient receiving the second identified treatment. This estimate is called the ratio of rate ratios. By using the weighted data in performing the regression, the estimate of the ratio of rate ratios is doubly robust.

FIG. 1 is a schematic block diagram of a medical treatment hypothesis testing system constructed according to an embodiment of the present invention for comparing the effects of medical treatments based on retrospective observational data. As shown in FIG. 1, the system includes an operator terminal 102 for operator access, over a network 104 to an analysis machine 106 and a patient data machine 108. The patient data machine 108 provides access for an operator to a patient record database 110 and an auxiliary database 112. The analysis machine 106 may include a processor; and memory operably coupled to the processor storing programming instructions and other data therein, the processor being operable to execute program instructions, a network connection to permit the analysis machine to receive input from the patient data machine 108 and other sources and to output results to the operator terminal 102 and other destinations.

The patient record database 110 stores retrospective observational patient data from patient medical records such as patient identification number, date of birth, gender, the results of medical tests, observational data recorded by healthcare providers, and information provided by patients to healthcare providers. The auxiliary database 112 stores information related to the operation of the system such as information about what records have been retrieved from the patient record database and when they were retrieved, system data formatting rules, and other data pertinent to the analysis of the patient record data.

The analysis machine 106 provides access for an operator to several software components including a data selection component 114, a propensity scoring component 116, a regression modeling component 118 and a data output component 120. These software components may take the form of computer instructions stored in computer memory and executed by a computer processor. The data selection component 114 presents the operator with an interface for the selection of relevant patient data and attributes from the patient record database 110, and retrieves and formats the selected data for use in the propensity scoring component 116.

The propensity scoring component 116 determines and assigns a propensity score to each patient record that represents that patient's likelihood of being in the group of patients receiving the identified treatment (“in the treatment group”). The propensity scoring component 116 further applies the propensity score to the patient data to weight the data of the patient records such that the weighted data for the group of patients not in the treatment group closely resembles the non-weighted data for the group of patients in the treatment group.

The regression modeling component 118 provides an interface for the operator to build and test a model for estimating the relative likelihood that a patient receiving the first identified treatment would experience an identified effect as compared to a patient receiving the second identified treatment (i.e. a patient not in the treatment group). The regression modeling component 118 receives weighted data weighted by the propensity scoring component 116.

The data output component 120 provides an interface to allow the operator to select the format and style for presenting analysis machine results. In some embodiments the data output component 120 includes a tool for selecting and formatting data produced by the analysis machine. In other embodiments the tools in the data output component 120 allow the operator to select and manipulate various visualization tools such as charts and graphs to assist interpretation and understanding of analysis machine results.

FIG. 2 is a high-level flow diagram of the steps exercised by the data selection component 114 of the medical treatment hypothesis testing system of FIG. 1 according to an embodiment of the present invention. As shown in FIG. 2, the data selection component 114 initially establishes a connection to the patient record database 202 to provide access to the patient data for the operator. According to one embodiment of the invention, the component utilizes an ODBC driver to connect the data selection component 114 to the patient record database 110 on the patient data machine 108 across the network 104. The component then provides a data selection interface to the operator 204 to allow the operator to select the appropriate patient records for the subsequent analysis on the analysis machine. According to one embodiment of the invention, the data selection interface is presented to the operator as an HTML webpage viewed through a web browser application. After the operator selects the desired data, the system receives this selection 206 and retrieves the selected records 208 from the patient record database 110. The system then formats the data 210 for use in the analysis machine by the propensity scoring component and stores the formatted data 212 in memory on the analysis machine 106.

FIG. 3 is a high-level flow diagram of the steps exercised by the propensity scoring component 116 of the medical treatment hypothesis testing system of FIG. 1 according to an embodiment of the present invention. As FIG. 3 depicts, the propensity scoring component 116 performs an initial data preparation step 302 described in greater detail with respect to FIG. 4. In this data preparation step 302 the analysis machine prepares the patient records in a table and identifies the initial “seeding” propensity score and residual associated with each patient record, creates columns for and evaluates each indicator function and the products of indicator functions for each patient characteristic described in the selected data. After the data preparation step 302, the analysis machine identifies the largest absolute correlation 304 between any indicator function column or product of indicator functions column and the residual column. For the following formulas, t_(i) represents the 0/1 indicator of patient i receiving the first identified treatment, p_(i) is the estimated probability that patient i received the first identified treatment, I_(j) is the j^(th) indicator function, and n is the number of patients. To determine the extent to which any two columns are correlated, in order to identify the largest absolute correlation, the following formula is employed: r_(j)=Sum_(i)((t_(i)−p_(i))(I_(ji)−mean(I_(j))))/((n−1)sd(t−p)sd(I_(j))). After identifying the column most correlated to the residual column, the analysis machine adjusts the propensity scores for all patients 306 based on the identified column. In adjusting the propensity scores the following formula is employed: p_(i)/(1−p_(i))=p_(i)/(1−p_(i))×exp(δ×sign(r_(j))×(2I_(ji)−1)), where δ is a tuning parameter set to a small number such as 0.001. Using the new propensity score that results from the adjustment process, the analysis machine then determines a weight to apply to each of the non-treated patient records 308. The weight is calculated from the propensity score using the following formula: w_(i)=p_(i)/(1−p_(i)).

The analysis machine compares the aggregate weighted data for patients not in the treatment group against the aggregate data for patients in the treatment group 309 to determine whether the two data sets are sufficiently similar (“optimally balanced”) 310. If the data sets are not optimally balanced, the analysis machine returns to the data table, that contains newly assigned propensity scores, assigned in the adjustment process 306 and the newly calculated residual (adjusted by the change of the propensity score). The analysis machine then repeats the process of identifying the largest absolute correlation 304 with the new residual, adjusting the patient propensity scores 306, determining new patient data weights based on the newer propensity scores 308 and comparing the weighted non-treatment data with the non-weighted treatment data 309 to determine whether the data sets are now optimally balances 310. This process repeats until the data sets are sufficiently similar, at which point the analysis machine formats the data for use in the regression step 312.

FIG. 4 is a flow diagram of the steps exercised by the Prepare Data in Table component of the propensity scoring component of FIG. 3 according to an embodiment of the present invention. The data preparation steps depicted in FIG. 4 interpret and arrange the selected patient data in a manner such that it can be interpreted by the subsequent step in the propensity scoring process involving correlation 304. The data preparation described in the foregoing steps detailed in FIG. 4 is also illustrated in exemplary tables in FIG. 5A-FIG. 5E. As FIG. 4 shows the system initially retrieves the patient case data 402 that was stored within the data selection component process 114 also depicted in detail in FIG. 2. An exemplary table containing patient case data for five example patients is illustrated in FIG. 5A. The illustration of FIG. 5A shows that each patient record at the initial stage may contain at least a patient ID number, a treatment indicator for identifying whether the patient was a member of the treatment group and data associated with categorical features such as race as well as numeric measurements such as blood pressure. Within FIG. 4, the next step in the process is the assignment of an initial “seeding” propensity score to each patient 404. In one embodiment of the invention, this seeding propensity score is calculated by the analysis machine dividing the number of patients within the treatment group by the total number of patient records. The next step in the data preparation process is the calculation of the residual for each patient record. In one embodiment of the present invention, the residual is calculated by subtracting the propensity score from the treatment indicator for each patient. Exemplary tables containing propensity score and residual data for each of the five example patients are illustrated in FIG. 5B and FIG. 5C respectively. Returning to FIG. 4, in one embodiment, the analysis machine then defines and evaluates indicator functions for each patient characteristic described in the data. The definition of said indicator functions, and creation of associated columns, is illustrated in the exemplary table of 5D that shows a new column for each condition that must be defined to allow the analysis machine to describe in boolean terms both categorical features and numeric measurements. In another embodiment the evaluation of all indicator functions is avoided by the use of published algorithmic shortcuts, in this case only selected indicator functions are defined and evaluated. The indicator functions are evaluated by determining whether the condition contained in the indicator function is true or false with respect to each patient record, with a 1 or 0 placed in the column to represent the same. The next step in the process of FIG. 4 is the definition and evaluation of the product of each indicator function with each other indicator function 410. As with the definition of the indicator functions, the system creates a new column for each condition and evaluates each condition placing a 0 or 1 in the relevant column to indicate whether the patient identified in that row meets the condition identified in the indicator function, in this case the product of one or more indicator functions, at the top of each column. The columns created by this step and the evaluation of patient records is illustrated in FIG. 5E.

FIG. 6 is an illustration of an exemplary table for the comparison of patient data for patients who received the treatment and weighted patient data for patients who did not receive the treatment according to an embodiment of the present invention. FIG. 6. illustrates the data considered in the comparison step 309 of FIG. 3 in determining whether the data has reached an optimal balance 310. The data in the exemplary table of FIG. 6 shows a sufficiently balanced data set after weighting to prepare the data for use in regression 312. The sufficiently balanced data set in the exemplary table of FIG. 6 has a maximum absolute difference less than or equal to 0.2 percent as between corresponding patient data of the aggregate weighted data for patients not in the treatment group and the aggregate data for patients in the treatment group.

FIG. 7 is an illustration of an exemplary patient case data table for use in the regression modeling 118 of FIG. 1. This illustration shows the data as it is formatted during the format data for regression step 312 of FIG. 3. In this example the original data from FIG. 5A is supplemented with data on the observed identified effects. The table of FIG. 7 shows two records for each patient, one for a “pre-treatment period” and one for a “post-treatment period.” The “period” column of FIG. 7 contains a 1 or 0 to indicate whether the record describes the “post-treatment” or “pre-treatment” period respectively. The table of FIG. 7 also contains an column/indicator “Y” that indicates whether the patient experienced the identified effect during the observation period, and a column/indicator “E” that describes the length of the observation period for each patient. The table of FIG. 7. also includes a column/indicator for propensity score weight, “W”, computed by the system 308 as depicted in FIG. 3 and described above.

The relative likelihood of a patient receiving the first treatment experiencing the identified effect when compared with a patient receiving the second treatment is calculated by estimating the ratio of rate ratios. The ratio of rate ratios (“RRR”) compares the rates of the identified event in the pre-treatment period to the post-treatment period for the treatment (first treatment) and non-treatment (second treatment) groups. The formula for calculating the RRR is: (rate(post,treatment)/rate(pre,treatment))/(rate(post,non-treatment)/rate(pre,non-treatment)). After calculation, the analysis machine outputs this estimated RRR to the operator terminal or another location where it may be used to evaluate the relative efficacy and risks of the first and second identified treatments. Decision makers in the medical field may utilize the estimated RRR output from the analysis machine to determine whether to approve or recommend particular treatments for patients and healthcare organizations.

To adjust for potential confounding the system obtains a doubly robust adjusted ratio of rate ratios by estimating a weighted Poisson regression model 118 as shown in FIG. 1. The estimates are derived from maximizing the equation L(b, δ)=Sum_(i)W_(i)(Y_(i)f_(i)−exp(f_(i))) where f_(i)=log(E_(i))+b₀+b₁T_(i)+b₂Period_(i)+b₃T_(i)*Period_(i)+δ′X. The resulting value of exp(b₃) is the doubly robust adjusted ratio of rate ratios.

After calculation, the analysis machine outputs the doubly robust adjusted ratio of rate ratios to the operator terminal for review by the operator, or to another location where it may be used to evaluate the relative efficacy and risks of the first and second identified treatments. Decision makers in the medical field may utilize this doubly robust adjusted ratio of rate ratios output from the analysis machine to determine whether to approve or recommend particular treatments for patients and healthcare organizations.

Although the present invention has been described through the use of exemplary embodiments, it will be appreciated by those of skill in the art that various modifications may be made to the described embodiments that fall within the scope and spirit of the invention as defined by the claims and their equivalents appended hereto. 

What is claimed is:
 1. An apparatus for evaluating effects of treatments comprising: a processor; and a memory operably coupled to the processor and having program instructions stored therein, the processor being configured to execute the program instructions, the program instructions comprising: receiving record data of a plurality of subjects; identifying relevant characteristics for evaluation; identifying a first treatment; identifying a second treatment; identifying an effect; identifying, from among the subjects, a first group of subjects that received the first treatment; identifying, from among the subjects, a second group of subjects who did not receive the first treatment; assigning a weight to the record data of each subject of the second group of subjects corresponding to a likelihood that the corresponding subject would have received the first treatment to generate weighted record data of the second group of subjects; determining whether the weighted record data of the second group of subjects is balanced with the record data of the first group of subjects; adjusting the weight assigned to the record data of each subject of the second group of subjects to balance the weighted record data of the second group of subjects with the record data of the first group of subjects; determining an estimated relative likelihood, using the weighted record data of the second group of subjects to perform a regression, that the first treatment, as compared to the second treatment, will result in the effect; and outputting the estimated relative likelihood.
 2. The apparatus of claim 1, wherein determining the estimated relative likelihood that the first treatment, as compared to the second treatment, will result in the effect further comprises using the record data of the first group of subjects, which is unweighted, in the regression.
 3. The apparatus of claim 1, wherein the likelihood that the corresponding subject would have received the first treatment is based on the relevant characteristics of the record data of the corresponding subject.
 4. The apparatus of claim 1, wherein assigning the weight to the record data of each subject of the second group of subjects comprises setting a propensity score and a residual value for the record data of each of the subjects.
 5. The apparatus of claim 4, wherein the record data comprises columns of indicator functions based on the relevant characteristics, each of the indicator functions being associated with a corresponding one of the subjects, wherein the residual value for the record data of each of the subjects is arranged in a residual value column, each place in the column corresponding to one of the subjects, and wherein assigning the weight to the record data of each subject of the second group of subjects further comprises: identifying which column of the indicator functions is most correlated to the residual value column; adjusting the propensity score for the record data of each of the subjects based on the identified column; and determining the weight to assign to the record data of each subject of the second group of subjects based on the propensity score.
 6. The apparatus of claim 5, wherein when identifying which column of the indicator functions is most correlated to the residual value column, the following equation is used: r _(j)=Sum_(i)((t _(i) −p _(i))(I _(ji)−mean(I _(j))))/((n−1)sd(t−p)sd(I _(j))).
 7. The apparatus of claim 5, wherein when adjusting the propensity score of the record data of each of the subjects based on the identified column, the following equation is used: p _(i)/(1−p _(i))=p _(i)/(1−p _(i))exp(δsign(r _(j))(2I _(ji)−1)).
 8. The apparatus of claim 5, wherein when determining the weight to assign to the record data of each subject of the second group of subjects based on the propensity score, the following formula is used: w _(i) =p _(i)/(1−p _(i)).
 9. The apparatus of claim 5, wherein the adjusting the weight assigned to the record data of each subject of the second group of subjects to balance the weighted record data of the second group of subjects with the record data of the first group of subjects comprises: adjusting the residual value of the record data of each of the subjects based on the propensity score; identifying which column of the indicator functions is most correlated to the residual value column; adjusting the propensity score for the record data of each of the subjects based on the identified column; and determining the weight to assign to the record data of each subject of the second group of subjects based on the propensity score.
 10. The apparatus of claim 4, further comprising setting a treatment identifier in the record data of each of the subjects, wherein the treatment identifier of the record data of each of the first group of subjects is set to 1, wherein the treatment identifier of the record of each of the second group of subjects is set to 0, wherein the propensity score of the record data of each of the subjects is initially set to be equal to the number of subjects in the first group of subjects divided by the total number of subjects in the record, and wherein the residual value of the record data of each of the subjects is initially set by subtracting the corresponding propensity score from the corresponding treatment identifier.
 11. The apparatus of claim 1, wherein the determining whether the weighted record data of the second group of subjects is balanced with the record data of the first group of subjects comprises: aggregating the record data of the first group of subjects to generate aggregate first data; aggregating the weighted record data of the second group of subjects to generate aggregate second data; and comparing the aggregate first data to the aggregate second data.
 12. The apparatus of claim 11, wherein the determining whether the weighted record data of the second group of subjects is balanced with the record data of the first group of subjects further comprises: determining that the weighted record data of the second group of subjects is balanced with the record data of the first group of subjects when the aggregate first data is substantially similar to the aggregate second data.
 13. The apparatus of claim 12, wherein the aggregate first data is substantially similar to the aggregate second data when a maximum absolute difference between corresponding ones of the relevant characteristics of the first aggregate data and the second aggregate data is less than or equal to 0.2 percent.
 14. The apparatus of claim 1, wherein the estimated relative likelihood that the first treatment, as compared to the second treatment, will result in the effect is determined by estimating a ratio of rate ratios.
 15. The apparatus of claim 14, further comprising adjusting the estimated relative likelihood that the first treatment, as compared to the second treatment, will result in the effect is determined by estimating a weighted Poisson regression model.
 16. A system for evaluating effects of treatments comprising: at least one record database coupled to a network, the at least one record database configured to store record data of a plurality of subjects; and a server coupled to the network, wherein the server is configured to: access, through the network, the record data of the plurality of subjects stored in the at least one record database; identify relevant characteristics for evaluation; identify a first treatment; identify a second treatment; identify an effect; identify, from among the subjects, a first group of subjects that received the first treatment; identify, from among the subjects, a second group of subjects who did not receive the first treatment; assign a weight to the record data of each subject of the second group of subjects corresponding to a likelihood that the corresponding subject would have received the first treatment to generate weighted record data of the second group of subjects; determine whether the weighted record data of the second group of subjects is balanced with the record data of the first group of subjects; adjust the weight assigned to the record data of each subject of the second group of subjects to balance the weighted record data of the second group of subjects with the record data of the first group of subjects; determine an estimated relative likelihood, using the weighted record data of the second group of subjects to perform a regression, that the first treatment, as compared to the second treatment, will result in the effect; and output the estimated relative likelihood.
 17. The system of claim 1, wherein determining the estimated relative likelihood that the first treatment, as compared to the second treatment, will result in the effect further comprises using the record data of the first group of subjects, which is unweighted, in the regression.
 18. The system of claim 17, wherein the estimated relative likelihood that the first treatment, as compared to the second treatment, will result in the effect is determined by estimating a ratio of rate ratios.
 19. An apparatus for evaluating effects of treatments comprising: a processor; and a memory operably coupled to the processor and having program instructions stored therein, the processor being configured to execute the program instructions, the program instructions comprising: receiving record data of a plurality of subjects; receiving relevant characteristics for evaluation; receiving an identification of a first treatment, a second treatment, and an effect; identifying, from among the subjects, a first group of subjects that experienced the first treatment; identifying, from among the subjects, a second group of subjects who did not experience the first treatment; assigning a weight to the record data of each subject of the second group of subjects according to a propensity score and a residual value of the record data of each of the subjects to generate weighted record data of the second group of subjects; determining whether the weighted record data of the second group of subjects is balanced with the record data of the first group of subjects; adjusting the weight assigned to the record data of each subject of the second group of subjects to balance the weighted record data of the second group of subjects with the record data of the first group of subjects; determining an estimated relative likelihood, using the weighted record data of the second group of subjects and unweighted data of the first group of subjects to perform a regression, that the first treatment, as compared to the second treatment, will result in the effect; and outputting the estimated relative likelihood. 